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We find new sum rules between direct CP asymmetries in D meson decays with coefficients that 
can be determined from a global fit to branching ratio data. Our sum rules eliminate the penguin 
topologies P and PA, which cannot be determined from branching ratios. In this way we can make 
predictions about direct CP asymmetries in the standard model without ad hoc assumptions on the 
sizes of penguin diagrams. We consistently include first-order SU(3 )f breaking in the topological 
amplitudes extracted from the branching ratios. By confronting our sum rules with future precise 
data from LHCb and Belle II one will identify or constrain new-physics contributions to P or 
PA. The first sum rule correlates the CP asymmetries ajf'p in D° —> K + K~, D° —> Tv + n~, and 
D° 7r°7r°. We study the region of the acp(D° —> TT + n~)-acp(D° —> n 0 n°) plane allowed by 
current data and find that our sum rule excludes more than half of the allowed region at 95% C.L. 

Our second sum rule correlates the direct CP asymmetries in D + —b K°I\ + , Df —> I\ 0 tt + , and 
D+ K+n°. 


INTRODUCTION 


Decays of charmed mesons are currently the only way 
to probe flavor violation in the up-quark sector. A major 
goal of experimental charm physics is the discovery of 
CP violation in nonleptonic charm decays (see, e.g., P- 
H)- To this end, it is promising to study singly Cabibbo- 
suppressed (SCS) decays d whose direct CP asymmetries 
may be large enough to be detected in the near future 


dir (d) = |A SCS (d)|MA SCS (d )| 2 

J |M scs (d)| 2 + |M scs (d)| 2 ' 


( 1 ) 


Here, A scs (d) is a D° decay amplitude and A scs (d) is 
the amplitude of the CP-conjugate D° decay. We write 
A scs (d) = A SGS (d) + A s b GS (d) with 


in terms of the strong phase 6(d) = arg [_4 s d(d)]. The 
smallness of ImA& renders aQp(d) highly sensitive to 
physics beyond the standard model (SM). To establish 
a “smoking gun” signal of new physics one needs reliable 
SM predictions for aQ T P (d) or at least robust theoreti¬ 
cal upper bounds on \a,£p(d) which cannot be exceeded 
within the SM. The difficulty of such theory predictions 
can be witnessed from Aa^lp = a^p(II 0 —> K + K~) — 
clqp(D 0 —> 7r + 7r _ ): estimates vary between 0(0.01%) [5] , 

0(0.1%) [Mini, -0.25% PH and-0.4% [J5], not 

excluding an enhanced SM value between ~ —0.6% and 
~ -0.8% [ME!]. There are claims that CP-violating ef¬ 
fects in charm physics can be 0(1%) ‘18]. The situation 
is not any better in CP violation induced by D—D mix¬ 
ing [15, . The key problem is our lack of knowledge of the 
penguin amplitude entering Ab(d) in Eq. (§ muni]. 


A s s S s (d) = X sd A sd (d ), Af cs (d) = - ^ A b (d), 


and the shorthand notation X q = V* * * § q V uq and A s d = 
(A s — Ad)/2 for the elements of the Cabibbo-Kobayashi- 
Maskawa matrix involved. The contribution of A b (d) to 
branching ratios B(d) is negligible, which essentially en¬ 
tails B(d) oc \ A s d(d)\ 2 . 

Adopting the Particle Data Group convention (with 
A s > 0) one may safely neglect subleading powers of 
Afe/Asd ~ 10- 3 . Eq. § reads 


a%p(d) 


ImA b lm[e~ iS ^A b (d)] 


( 2 ) 


All theoretical analyses of D —> PP' decays, where 
P,P' denote pseudoscalar mesons, rely on the approxi¬ 
mate SU (3 )p symmetry of the strong interaction [7J flPb 
HI (ME]- In analyses of branching ratios one can in¬ 
clude first-order SU(3 )p breaking [71 HUHT71[?fl 1551151111521 
[331 [35U37U381SQUH] . An intuitive way to exploit SU(3 )f 
relations involves topological amplitudes (pioneered in 
Refs. [551 ES Sl|) which characterize the flavor flow in 
terms of tree (T), color-suppressed tree (0), exchange 
(E), annihilation (A), penguin ( Pd, s ,b ), and penguin an¬ 
nihilation ( PAd, s ,b ) diagrams. This method has been ex¬ 
tended to include linear SU (3) f breaking in applications 
to B [45| and D 01] decays. The first-order SU(3 )f- 
breaking corrections are parametrized by C^A^,.. 
with i = 1,2,3 labeling which d-quark line is replaced by 
an s line. 
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With our inability to predict individual CP asymme¬ 
tries it is natural to study correlations among several 
asymmetries. There are two sum rules which hold in the 
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limit of exact SU(3 )f symmetry J7J: 161 [311: 

a d j P {D° -> K+K~) + a A S T P {D° -A tt+tt" ) = 0, (3) 
a^ P (D+ -A K°K+) + a d j P {Dt -A K°tt+ ) = 0 . (4) 

Analyses of branching ratios permit the determination of 
|Af^ s (d)| in Eq. ([ 2 ]) and, through global fits, also to con¬ 
strain the phase 5(d). The branching ratios of the decays 
entering Eqs. © and 0 exhibit sizable SU(3 )f breaking 
which limits the power of these SU(3)j?-limit sum rules 
to test the SM. In this Letter, we derive new sum rules 
which incorporate SU(3 )f breaking in A®^ s (d) to linear 
order. To this end, we use the result of our global fit to 
D -A PP' branching ratios in Ref. [IT] in two ways: On 
one hand, we extract |Af)p s (d)| and 5(d) for the decays 
of interest from the fit to find the numerical relation be¬ 
tween acp(d) and Ab(d) in Eq. 0. On the other hand, 
the same fit also returns the values of the topological 
amplitudes entering the fitted branching ratios. How¬ 
ever, the desired Ab(d) of individual decay modes involve 
additional topological amplitudes (new combinations of 
Pd.s.b and PAd, s ,b) which do not appear in the branch¬ 
ing ratios. Our sum rules are constructed in a way to 
eliminate these unknowns. Unlike Eqs. © and © , these 
new sum rules use the SU(3)p limit for the eliminated 
Pd, s ,b , PAd, s ,b hr Ab(d) only, while consistently including 
SU(3 )f breaking in the SM-dominated quantities T, A, 
C, E. Since no relations among direct CP asymmetries 
hold to first order in SU(3 )f breaking [44], this is the 
best which can be achieved. 


with cf specified in Table [I] Since e~ zS ^ A s d(d) = 
|A s d(d)| is real, c d d A s d(d) does not contribute to a,Qp(d) 
in Eq. 0. 

Our strategy involves three steps: In step 1, we de¬ 
termine all quantities entering a cp(d) in Eq- 0 except 
for P and PA from a global fit to branching ratio data 
as described in Ref. [3Tj . These fitted quantities are the 
topological amplitudes A^\ C, C^\ and Pbreak entering 
Eq.0 as well as \A s d(d)\ and the phases 5(d). Note that 
|A s d(d)| is trivially found from B(d) oc \A s d(d)\ 2 . As an 
important feature of charm physics, the large number of 
different branching fractions gives useful information not 
only on the magnitudes of the topological amplitudes, 
but also on their phases (up to an overall unphysical 
phase). By plugging the results back into A s d(d), we 
find 5(d) (up to discrete ambiguities). 

In step 2 we eliminate all hadronic parameters but P 
and PA from a^)p(d). To this end, we define 


S(d) = 2 ie l ^ d) 


' \AT(d)\a d j P (d) 

Im Af, 


ImA(d) , 


(9) 


with X(d) given in Table [IT} S(d) is calculated solely from 
experimental input, as all ingredients of Eq. 0 are found 
from the global fit to branching ratios. To relate S(d) to 
theoretical quantities, we use Eq. © 


S(d) 


2*e i5(d) Im 


e~^ d) A b (d) - X(d) 


( 10 ) 


CP ASYMMETRY SUM RULES 


For our example above the subtraction term -ImA(rf) 
removes Pbreak and Eq. (10) gives 


Our analysis starts with the decomposition of A s d(d) 
and Ab(d) in terms of topological amplitudes. As an 
example consider 

Asd(D° -A 7T+7T-) = -T - E + Pbreak (5) 

Ab(D° -a tt+tt") = T + E + P + PA (6) 

with the SU(3)F-breaking penguin topology Pbreak = 
P s — Pd and s, d denoting the quark flavor in the pen¬ 
guin loop. Ab involves the new penguin topologies P = 
Pd + P s — %Pb and PA = PAd + PA S — 2PA&, which do 
not appear in A s d(d) of any decay d and, consequently, 
cannot be constrained by data on branching ratios. Next 
we combine Eqs. 0 and © to eliminate the numerically 
large parameters T and E from Ab 

Ab(D° —► TT+TT") 

= —A s d(D° —> 7T+7T ) + Pbreak + P + PA . (7) 

Trading T and E for A s d is possible for all decays d 

considered in this Letter: With 7) = C, Pbreak, •, we 
write 


S(D° — > 7T + 7T ) 

= 2ie i5{D °^ 7I+ *~ 


^Im 


\P + PA) . (11) 


The right column of Table [TT] shows S(d) in terms of P, 
PA, their complex conjugates P*, PA* , and 5(d) (which 
is determined from the fit) for all decays. (To relate 
Eq. @ to the entry in Table |IlJ use 2*Imz = z — z* for 
Z = e -i<5(£>°^7r+7T-)(p + PA).) 

In step 3 of our analysis we construct two sum rules 
eliminating P, PA, P*, and PA*, each of which con¬ 
nects three a^ T P (d). Table |ll| reveals that the combined 
information from S(D° —► tt^t:~) and S(D° —► 7r°7r°) de¬ 
termines the complex quantity P + PA. Any additional 
CP asymmetry depending on P + PA will then probe the 
standard model; i.e., with three CP asymmetries we can 
construct the desired sum rule eliminating P + PA: 

CP Asymmetry Sum Rule 1: 

S(D° K+K~) - S(D° ->• tt+tt-) 

e 2iS(D°^K+K-) _ e 2i5(D°-+Tr+ir-) 

S(D° -A K+K~) + V2S(D° -A 7r°7T 0 ) 


A b (d) = c d sd A sd (d)+ Y] ° d iTi , 


( 8 ) 


e 2i8(D°^K+K-) _ e 2i8(D°->TT 0 TT°) 


0 . ( 12 ) 
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Decay amplitude A(d) A , 

,d(d) A' 1} 

4 1} 

4 15 

C + A cf } 

-Pbreak P +2 A P~\~P A 

A(D° -A K+K-) 

1 

0 

0 

0 

0 

0 

-l 

0 

i 

A(D° —> 7T + 7T“) 

-1 

0 

0 

0 

0 

0 

l 

0 

l 

A(D° -A 7r°7r°) 

-1 

0 

0 

0 

0 

0 

1 

a/2 

0 

1 

a/2 

A{D+ -A- K°K+) 

1 

0 

0 

2 

0 

0 

-l 

1 

0 

A{D+ -A K° 7T+) 

-1 

2 

2 

0 

0 

0 

l 

1 

0 

A(D+ -A K+ 7T°) 

1 

0 

0 

0 

V2 

y/2 

1 

a/2 

1 

a/2 

0 


TABLE I. The coefficients cf of the topological decomposition of Ab{d) in Eq. (jsj) . 4f 3 and cf' are first-order SU(3 )f- 
breaking corrections to A and C, respectively, as defined in 1411 . 


Decay d 

X(d) 

S(d) 

D° -a K+K~ 

e- iS ^ U ^ K +K ~H- P ^) 

(P + PA)-e 2iS ( DU ^ K+K )(P + PA)* 

D° -A n + n~ 

(P break ) 

(P + PA) - e ™( DU ^ + ”~) (p + PA)* 

7~)0 . _0 0 

U — >■ 7T 7T 

e- ,j ( D ^ ,0 ’°) (-fjPbreak) 

W (_p _ PA) _ e 2iS{n°^n°n°) _ pA y 

D+ -a K°K+ 

e -H D+ ^ KU K + ) (-2 A taC + ^ R0K+ + 25a ~ Pbreak) 

p g 2 iS(D+^K u K + ) p* 

D+ -A K° 7T+ 

e -« (D +^^+) (2A^ k0 ^ + 25 a + P break ) 

P _ g2i«(D+->if°vr+)p* 

D+ -A K + 7T 0 

e -iS{n+-*K+*“) + ^2 C m + ^p break ) 

^(-P)-e 2iS{Dt " K+nU )^(- p r 


TABLE II. Definitions of X'(d) and results for < 5(d) as used and defined in Eq. |9]) in case of the SU(3 )f fit including 1/N C 
counting. For the sign conventions of A fac (d) see Ref. |41| . Note that A 1 p%^ t ^a K + = 0 by isospin symmetry. 


In this last step SU(3 )f breaking in P, PA is neglected. 
The second sum rule, which correlates three CP asym¬ 
metries depending on P, is 

CP Asymmetry Sum Rule 2: 

S(D+ -A K°K+) - S{D+ -A K°n+) 

e 2i5(D+^K°K+) _ e 2i6(D+^K°-K+) 

S(D+ -A K°K+) + y/2S(D+ -A K+n°) _ 

e 2iS(D+^K°K+) _ e 2i5(D+^K+Tr°) ' ^ ' 


By inserting the expressions in the right column of Ta¬ 
ble [IT] one readily verifies Eqs. (12) and (13). If some of 
the phases in the denominators of the sum rules are equal 
(covering the SU(3 )f limit as a special case) one finds: 
for e 2lS ( D ~* K+K ) = e 2lS ( D ->-n-+7r ) Eq. (121 collapses 
to 


S(D° -a K + K~) — S(D° -a 7r + 7r _ ) = 0, (14) 

while for e 2lS ( D °^ K+K ) = e 2 iS(D 0 ^>A\ 0 ) sum rl j e 
becomes 


S(D° -A K+K-) + V2S{D° 


mo 


) = 0. (15) 


If all three phase factors are equal, Eqs. (14) and (15) 
hold simultaneously. The special cases of sum rule 2 are 
obtained from those in Eqs. (14) and (15) by obvious 
replacements. 

When linking sum rule 2 to experimental quantities, 
we use a^p{D+ -a K S K+) = a$ r P (D+ -4 K°K+) and 


a cp{D+ —► Kstt + ) = a,Qp(D+ —» A'°7t + ), meaning that, 
in our definition of a^p kaon CP violation is properly 
subtracted [45]. The two sum rules probe the SU(3)f 
limit in P and P + PA. If future experiments find de¬ 
viations of order 30%, one will ascribe those to SU(3)f~ 
breaking hadronic effects. The smallness of Im Ab makes 
the sum rules highly sensitive to new physics, which may 
well violate the sum rules at a far higher level. 


SM PREDICTION OF CP ASYMMETRIES 


We combine the sum rules Eqs. (12) and (13) with 
the branching ratio fit presented in Ref. [41]: for each 
point in the parameter space complying with all mea¬ 
sured branching fractions (and the strong phase 5fCn), 
we determine S(d) for the decays entering the sum rules. 
In the same step Eqs. , and © are used to pre¬ 

dict one CP asymmetry in terms of the other two. In our 
fit, we demand \(C+6 A )/T iac \,\(E+6 A )/T iac \ < 1.3 to 


enforce proper l/N c counting. (T fac and A iac are the 
factorized tree and annihilation amplitudes and S A = 
A(D+ -A K°n+) - A fac (£>+ IC°tt + ) = 0(1/N 2 ). The 


very conservative choice <1.3 accounts for a large Wil¬ 
son coefficient in E, C offsetting the suppression factor 
of 1/-/V C ~ 0.3 [41] .) Apart from the fit with current data 
(see Table X of Ref. [H]), we also consider a hypothet¬ 
ical future scenario with improved branching ratios by 
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Observable Measurement References 


A acp 

Sagj, 

a d c %(D° -4 KsKs) 
acp(D° -4 W 0 ) 
a^(D+ -4 7r°7r+) 
agj»(D+ -4 KsK + ) 
a d j P (D+ -4 A s7 r+) 
agf,(ZP+ -4 A ”*"7r°) 


—0.00253 ± 0.00104 
—0.0011 ± 0.0026 
-0.23 ±0.19 
-0.0004 ±0.0064 
±0.029 ± 0.029 
±0.0011 ±0.0017 
±0.006 ±0.005 
±0.266 ±0.228 


HMS6H52 
f m im SZ1 Em E2] 
ESI 
+ mi53] 

EH 

f [5111551 l55H5n| 

EH 


TABLE III. Current data on SCS charm CP asymmetries 
with subtracted indirect CP violation from kaon and charm 
mixing [45l 161] , see Appendix A of Ref. [16] . We use the no¬ 
tation E a'g'p = aop(D° -4 K+K~)+ a^ P {D° -4 n+Tv~). No 
correlations between CP asymmetries are taken into account 
in the fits. ^Our average. Table adapted from Ref. [62j . 


branching ratios in our future scenario do not improve 
the predictions for CP asymmetries. Note that, with 
current data, the predicted ranges barely depend on the 
additional input from 1/N C counting. 

In principle, one can obtain the quantities P and P + 
PA, which we eliminate through our sum rules, from 
the individual CP asymmetries. If future data challenge 
our sum rules at a level which cannot be explained with 
SU(3) f breaking in P and P ± PA, this will point to new 
physics which couples differently to s and d quarks. (For 
an SU(3 )p analysis of such models see, e.g., Ref. [IS].) 
However, SU(3)F-symmetric new physics in P and P + 
PA vanishes from the sum rules. (A similar situation can 
be found in the isospin sum rules of Ref. [63] which are 
insensitive to new physics in A7 = 1/2 amplitudes.) 


scaling their errors with a factor l/-\/50. To illustrate 
the impact of the 1 / N c counting for the SM predictions, 
we perform an additional fit without 1/N C input. This 
plain SU(3 )f fit relies on the topological parametriza- 
tion of Table III in Ref. [41] with the SU(3 )f counting 
described in Sec. IIIB of Ref. EH- The redundancy of the 
four SU(3) p-limit topologies [4T] is removed by absorbing 
A into T, C, and E. We further demand |A^/T| < 50% 
to respect the SU(3 )f counting. 

The experimental values of the CP asymmetries in¬ 
cluded in the fit are summarized in Table m Our global 
fit results are shown in Fig. [l] The % 2 of the global min¬ 
ima range from 0.0 to 2.0 in the considered scenarios, i.e., 
with or without 1/N C counting and with current or fu¬ 
ture data, indicating an excellent fit. The sum rules have 
nontrivial implications for direct CP asymmetries, espe¬ 
cially when combined with the input from 1 /N c counting. 
With current data, we see the largest impact of the sum 
rules in Fig. |l(a)| roughly 47% of the 95% C.L. region 
allowed by the measurements of the CP asymmetries is 
excluded by our global fit result. 


Our results show that future improved measurements 
of branching ratios will play a key role to sharpen 
our predictions: drastic exa mples are the prediction of 
aQp(D+ -4 K + tt°) in Figs. 1(c) and 1(d) and the cor¬ 
relation of clqp(D° -4 7r°7r°) and TF tt+tt") in 
Fig - [ H aj| where one of the two overlapping ellipses van¬ 
ishes in our future-data scenario. In Fig. |l(b)[ a smaller 
portion of th e expe rimentally allowed region is excluded 
than in Fig. 1(a) because a^p^D 0 —>■ 


is mea¬ 


sured less precisely than Aa^p (see Table |Tu| ), rendering 
the sum rule less powerful in Fig. |l(b)| In general, the 
correlation of two CP asymmetries can be better pre¬ 
dicted once improved data for the third one appearing 
in Eqs. (121 and (13) becomes available. In the case of 
the fit without l/iV c counting the smaller errors of the 


CONCLUSIONS 


To find reliable SM predictions for charm CP asym¬ 
metries, we derive two sum rules which treat T, A, C, 
E, Tbreak correctly to linear order in SU(3 )f breaking 
and eliminate the penguin topologies P and P ± PA to 
leading order in SU(3 )f- Thus, we treat large tree-level 
parameters at subleading order to increase the precision 
in the extraction of loop-induced quantities sensitive to 
new physics. Unlike previously known SU(3)F-limit sum 
rules, our new sum rules correlate three CP asymmetries 
each. The interplay of the two sum rules probes both 
the quality of SU(3 )f for penguin topologies and new 
physics. Future branching ratio measurements play a key 
role in order to reduce the uncertainties of the consequent 
SM predictions for charm CP asymmetries. 
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(a) Fit including 1/N C counting in the topological amplitudes 
with current (future) branching ratio data in blue (green). The 
black (magenta) line delimits the 95% C.L. region found in a 
fit without using 1/N C counting in the topological amplitudes 
for current (future) branching ratio data. Note that the black 
and magenta curves lie on top of each other. 


(b) Color coding as in Fig. |l(a)] The dash-dotted line denotes 
the SU(3)_F-limit sum rule of Eq. § 



a$' p (Dt -> K+, r») 


(c) Prediction from current branching ratio data. The dashed, 
solid, and dash-dotted lines correspond to the la, 2a , and 30- 
intervals, respectively. The blue bars show our fit results for 
current branching ratio data. The corresponding result without 
1/N C counting (i.e. only using SU(3 )f) is shown in black. 



(d) Same as in Fig. (c), but for our future scenario with 
smaller errors of the branching ratios. The green (magenta) 
bars correspond to the analysis with (without) 1/N C counting. 



(e) Future scenario assuming a^p(D^ K+tt°) = — 1% and 
branching ratios with current (future) errors including 1/N C 
counting in blue (green). The corresponding allowed area 
found in a fit without 1/N C counting is delimited by the black 
(magenta) line. Note that the black, magenta, and red curves 
lie on top of each other. The dash-dotted line denotes the 
SU(3)jr-limit sum rule of Eq. 01- 


FIG. 1. SM predictions for CP asymmetries obtained from our global fit. In Figures (a),(b),(e) the dashed (solid) red lines 
delimit the experimental 68% (95%) C.L. regions. The other dashed (solid) lines are the 68% (95%) C.L. regions of the 
respective fit scenarios explained in the captions of the subfigures. The generic error of order ~ 30% from SU(3 )f breaking in 
P and P + PA is not shown. The experimental error ellipses are obtained by scaling the errors quoted in Table m by factors 
of y/2.28 and V5.99 in order to obtain the two-dime nsiona l 68% and 95% C.L. regions from the corresponding one-dimensional 
ranges. The experimental error ellipse shown in Fig. 1(b) is calculated from the corresponding one for A a^p and SaJ-jp. 
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